Observation of Superfiuid Flow in a Bose-Einstein Condensed Gas 
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We have studied the hydrodynamic flow in a Bose-Einstein condensate stirred by a macroscopic 
object, a blue detuned laser beam, using nondestructive in situ phase contrast imaging. A critical 
velocity for the onset of a pressure gradient has been observed, and shown to be density dependent. 
The technique has been compared to a calorimetric method used previously to measure the heating 
induced by the motion of the laser beam. 
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Beginning with the London conjecture [pj, Bose- 
Einstein condensation has been considered crucial for the 
understanding of superfluidity. Since then, the weakly in- 
teracting Bose gas has served as an idealized model for 
a superfiuid [H . It has a phonon-like energy-momentum 
dispersion relation that does not allow for the generation 
of elementary excitations below a critical velocity, thus 
implying dissipationless flow at lower velocities. The on- 
set of dissipation has been treated in the framework of 
the nonlinear Schrodinger equation ||-f7| and shows an 
intriguing richness. Experiments on liquid helium could 
not test these theories, because superfluidity and dissi- 
pation are even more complex in this system due to the 
presence of strong interactions and surface effects M . 

The creation of Bose-Einstein condensates in dilute 
gases has dramatically changed this situation, allowing 
for quantitative tests of microscopic theories using the 
tools and precision of atomic physics experiments 
A number of recent experiments have examined phe- 
nomenological features of superfluidity in gaseous Bose- 
Einstein condensates. These include the observation of 
vortices [ pO[[Tl] ] , a non-classical moment of inertia p2| , 
and suppression of collisions from microscopic impurities 
|jl3f . In previous work we found evidence for a critical ve- 
locity in a stirred condensate Q . In this Letter we study 
the onset of dissipation with higher sensitivity using re- 
peated in situ non-destructive imaging of the condensate. 
These images show the distortion of the density distribu- 
tion around the moving object, thus directly probing the 
dynamics of the flow field that has been recently treated 
with different models [fi"5|-p0|] . The experimental setup 
was similar to the one used in our previous work @. Im- 
provements in the evaporation strategy and decompres- 
sion techniques allowed us to produce pure condensates 
with up to 5xI0 7 sodium atoms, with densities rang- 
ing from 8.4xl0 13 to 3.5xl0 14 cm~ 3 , corresponding to 
chemical potentials from 60 to 250 nK. We determined 
the Thomas-Fermi radius R z along the axial direction 
through in situ phase contrast imaging. The sound ve- 
locity at the center of the condensate was then evaluated 
through the relationship c s = 2ttv z R z / 'V%, with v z =20.l 



Hz being the axial trapping frequency. The macroscopic 
moving object was a 514 nm laser beam blue-detuned 
with respect to the sodium transitions, thereby creating 
an effective repulsive potential for the atoms. The beam 
is focused on the center of the condensate to a Gaus- 
sian 1/e 2 diameter of 2w = lOyum. It is scanned with an 
acousto-optic deflector driven by a triangular wave along 
the axial direction of confinement of the magnetic trap. 
This results in a motion of the laser beam at constant 
speed v = 2a f, where / is the frequency and a is the am- 
plitude of the scan. To avoid effects due to the inhomo- 
geneous density distribution in the condensate, the laser 
beam was scanned at various frequencies with a constant 
amplitude approximately equal to the Thomas-Fermi ra- 
dius R z of the condensate along the axial direction. 
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FIG. 1. Hydrodynamic flow in a Bose-Einstein condensate. 
Eleven phase contrast images of a condensate stirred for one 
cycle are taken in situ for various velocities of the laser beam. 
From top to bottom: 0.4 mm/s, 7.0 mm/s, and 11.2 mm/s. 
The exposure time is 300 [is per frame, and the frame rates are 
33 Hz, 526 Hz, and 909 Hz, respectively. The maximum den- 
sity and the corresponding sound velocity without the laser 
are 2 x 10 14 cm -3 and 7 mm/s, respectively. The position of 
the laser beam is seen as a depletion of the condensate around 
its center. The vertical size of each image is ~ 200 /jm. 

The frequency range was limited to 200 Hz. Above 
this frequency the atoms experience an effective time- 
averaged stationary potential. The condensate was im- 
aged using the phase contrast technique fl9pl|] with probe 
light red-detuned from the F=l to F'=2 transition. Dif- 
ferent detunings (0.4 and 1.7 GHz) allowed us to maintain 
a large linear dynamic range when imaging condensates 
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at different densities. Fig. |l| shows multiple phase con- 
trast images of the hydrodynamic evolution of the density 
distribution of a single condensate during the stirring. 
For low velocities there was little effect on the conden- 
sate. However, at higher velocities the density distribu- 
tion became asymmetrical and the condensate piled up 
in front of the moving beam, resulting in turbulent flow 
around the object. 

The main result of this Letter is based on quantitative 
analysis of images like those in Fig. |l|. The analysis of 
the observed flow pattern shows a critical velocity for the 
onset of a drag force between the macroscopic laser beam 
and the condensate. 
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FIG. 2. Pressure difference across a laser beam moving 
through a condensate. On the left side in situ phase con- 
trast images of the condensate are shown, strobed at each 
stirring half period: beam at rest (top); beam moving to the 
left (middle) and to the right (bottom). The profiles on the 
right are horizontal cuts through the center of the images. 
The stirring velocity and the maximum sound velocity were 
3.0 mm/s and 6.5 mm/s respectively. 

For a weakly interacting Bose-condensed gas at den- 
sity n(r) and chemical potential /i, pressure is identical 
to the mean-field energy density P = /x(r)n(r)/2 [||. A 
drag force arises due to the pressure difference across 
the moving object. The chemical potential is given by 
n(r,t) = gn(r,t), where g = 4ir?i 2 a/M is the strength of 
two-body interactions for atoms with mass M and scat- 
tering length a. Therefore the magnitude of the drag 
force F is given by 

F ~ gSnAn = SfiAfj,/g (1) 

where An and A/x are the differences in density and 
chemical potential across the stirring object, and S the 
surface area the object presents to the condensate. 

If the laser beam is stationary, or moves slowly enough 
to preserve a possible superfluid state for the condensate, 
there will be no gradient in chemical potential across the 
laser focus, and therefore zero force according to Eq. ([!]). 
A drag force between the moving beam and the conden- 
sate is indicated by an instantaneous density distribution 
n(r, t) that is distorted asymmetrically with respect to 



the laser beam. Fig. g shows this asymmetry as the for- 
mation of a bow wave in front of and a stern wave behind 
the moving laser beam. 
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FIG. 3. Drag force on a Bose-Einstein condensate, derived 
as an asymmetry in the density profile (top). Beyond 4 mm/s 
the signal saturated due to the formation of a high-density 
bow wave, introducing nonlinearity into the phase contrast 
signal for phase shifts exceeding n. The bars represent sta- 
tistical errors. Scanning the laser at the axial trapping fre- 
quency, corresponding to a stirring velocity of 1.5 mm/s for 
an amplitude of 37 /im, shows synchronous oscillations of the 
condensate (bottom) and reduced heating. 

To disentangle asymmetries due to the moving beam 
from the intrinsic inhomogeneity of the density profile we 
strobed the images at half the stirring period and always 
viewed the laser beam in the center of the condensate. 
We defined the asymmetry A as the relative difference 
between the peak column densities in front (ri/) and be- 
hind (h b ) the laser beam A = 2(h f - h b )/(h f + h b ). We 
averaged A over 9 successive images of a single conden- 
sate. 

The average asymmetry versus velocity is shown in Fig. 
|. Below - 0.6 mm/s the asymmetry was zero within the 
statistical error indicating a suppression of dissipation in 
the superfluid regime. The asymmetry increased for ve- 
locities greater than 0.6 mm/s. The dip near 1.5 mm/s 
corresponds to a stirring frequency of ~ 20 Hz for the 
chosen amplitude. This is close to the axial frequency of 
the magnetic trap. Near this frequency the laser beam 
synchronously excited the axial dipole motion of the con- 
densate (see Fig. |[), thus reducing the relative velocity 
of the stirrer and the condensate and therefore the drag. 
We repeated the measurements at different stirring am- 
plitudes to disentangle velocity and frequency dependent 
effects. This confirmed that the observed dip always oc- 
curred near 20 Hz, whereas the other features of the ob- 
served asymmetries depended only on velocity. 
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FIG. 4. Density dependence of the critical velocity. The 
onset of the drag force is shown for two different condensate 
densities, corresponding to maximum sound velocities of 4.8 
mm/s (•, left axis) and 7.0 mm/s (x, right axis). The stirring 
amplitudes are 29 /im and 58 fim, respectively. The two ver- 
tical axes are offset for clarity. The bars represent statistical 
errors. 

To study the superfluid regime at different atomic den- 
sities, we focused on the low velocity region (v < 2 mm/s) 
by choosing the stirring amplitudes in such a way that 
the frequencies were below the trap resonance. The den- 
sity was varied by adjusting the number of atoms for 
constant trapping parameters. The power in the laser 
beam was scaled to keep the ratio of the potential height 
of the laser to the chemical potential equal to 8 within 
10%, thus preserving the effective size of the laser beam. 

In Fig. |] we show measurements of the average asym- 
metry for two maximum densities no of 9 x 10 13 and 



1.9 x 10 i4 cm- 



Two features are evident. In each 



data set there is a threshold velocity v c below which 
the drag force is negligible, and this threshold increases 
at higher density. Above this critical velocity, the drag 
force increases monotonically, with a larger slope at 
low density. The data were fitted to a piecewise lin- 
ear function that reflects the two regimes of dissipation, 
one with zero asymmetry below v c and another, above 
v c , where the asymmetry is linear in v — v c . The fits 
yield v c = (0.11 ± 0.02)c s for the higher density and 
v c = (0.07 ± 0.02)c s for the lower density. This is in 
agreement with the expectation that the critical veloc- 
ity increases with the speed of sound and disagrees with 
models that predict only a dependence upon the diameter 
of the object E^jl§|. 



One can compare measurements of the asymmetry 
(proportional to the drag force F) with the power trans- 
ferred to the condensate by the moving beam, F ■ v, us- 
ing an improved implementation of the calorimetric tech- 
nique introduced in |]l4j| . For this, we stirred the conden- 
sate for times between 100 ms and 8 s, in order to pro- 
duce approximately the same final temperature. After 
the stirring beam was shut off, the cloud was allowed to 
equilibrate for 100 ms. The thermal fraction was deter- 



mined using ballistic expansion and absorption imaging 
P,H . We inferred the temperature and total energy us- 
ing the specific heat evaluated as in 1 22 1. A further im- 
provement was obtained by comparing successive images 
with and without stirring of the condensate, to subtract 
out all other sources of residual heating. As a result of 
the improved fitting and background subtraction proce- 
dure, we were able to detect changes of less than 10 nK in 
the energy of the gas, as we will describe in more detail 
in p3| |. A condensed fraction higher than 90% ensured 
that the laser beam primarily heated the condensate and 
not the thermal cloud. 
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FIG. 5. Calorimetry of a condensate. The energy transfer 
rate during stirring (•, left axis) was obtained from temper- 
ature measurements. The error bars reflect shot-to-shot vari- 
ations in the temperature. The results are compared to the 
energy transfer rate vA(v) obtained from the asymmetry data 
of Fig. |^ ( x , right axis) . 

Fig. H shows the energy transfer rate to the condensate 
versus the stirring velocity. Here, as in the phase-contrast 
measurements, the amplitude of the stirring was kept 
fixed while scanning the frequency. We see the onset of 
dissipation near the value obtained using the drag force 
method, v c ~ 0.5 mm/s = 0.1 c s . 

The calorimetric measurements can be compared with 
the drag force inferred from the asymmetric density dis- 
tribution. Using Eq.(|l|), the energy transfer rate per 
atom is written in terms of the asymmetry as 



dE 
~dt 



F ■ v fJ^onod 



N 



N 



vA(v) 



(2) 



Here we have used S = l z d for the effective surface area of 
the object, with d the diameter of the laser beam, l z the 
effective column depth, and nn — hqI z for the observed 
column density. 

The comparison between the calorimetric and the drag 
force measurements is also shown in Fig. |^, with an ad- 
justed scale factor for vA(v). The measurements of force 
(as an asymmetry in the instantaneous condensate wave- 
function) and heating (as an integrated increase in ther- 
mal fraction) are in remarkable agreement and demon- 
strates the consistency between our two methods of an- 
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alyzing dissipation. For the parameters of our experi- 
ment (d ~ 10/im, no = 1-3 x 10 14 cm -3 , l z — 66//m, 
N = 1.8- 10 7 ) the overall heating rate predicted by Eq.(||) 
is 4.5 times larger than that obtained from calorimetry. 
This discrepancy is not surprising, because Eq.(Q) does 
not account for dynamical pressure, the inhomogeneous 
density profile, and the 3D nature of the flow field. 

In one model ||, the drag force arises from the peri- 
odic emission of vortex lines at a rate that increases with 
velocity. The vortices reduce the pressure gradient across 
the object and the predicted heating rate is 0,0] 



dE 1 fj, v . . 



(3) 



where e VO rtex — ^irnfi 2 /Mla(d/£) is the vortex energy 
and £ = (8nna) is the healing length. Equating this 
to the heating rate in Eq.(Q), predicts that the slope of 
the asymmetry a = A(v)/(v — v c ) should scale as a ~ 
1/Cg. For the two densities in Fig. [|, the model predicts 
the ratio ai ow / cthigh = 2.1, while our measurement yields 
(1.9 ± 0.2). This, and the absolute value of the observed 
heating rate (see Ref . |b|] ) , are further evidence that the 
drag force is closely linked to vortex shedding. 

Our new measurements are consistent with our ear- 
lier experiment fl4|| where dissipation at subsonic veloc- 
ities was observed. However, the scatter in the observed 
heating rates in |lj] prevented a clear observation of 
the threshold, and the critical velocity was determined 
by linear extrapolation from high heating rates to be 
v c /c s = 0.25 ± 0.04. Other possible extrapolations tend 
to lower the value of v c , for instance a function cx v{v—v c ) 
as suggested by Eq.(§), yielded v c /c s = 0.17±0.09. With 
improved calorimetry we have taken data for parameters 
similar to |14| in the range 0.1 < v/c s < 0.25 and ob- 
served a small heating rate (~ 20 nK/s at 1 mm/s) that 
was previously indiscernible. 

Since our first report on critical velocities in a Bose 
Einstein condensate, several theoretical papers have 
added further insight p^-pcfl . The observed critical Mach 
number v c /c s ~ 0.1 is lower than the predictions for ho- 
mogeneous |rj,|5| and inhomogeneous || 2D systems. This 
discrepancy is most likely due to the 3D geometry of our 
experiment, where the laser beam pierced lower density 
regions of the condensate. It was noted in Ref. |{l8| that 
the critical velocity for phonon excitation is lowered by 
the inhomogeneous density distribution. Ref. jl6) shows 
that vortex stretching and half-ring vortices can lower the 
3D critical velocity below the 2D value. Jackson et al. 
fL5[ performed 3D simulations for a geometry similar to 
our experiment and obtained a critical velocity as low as 
0.13c s , in good agreement with our results. The relevant 
critical velocity in our experiment is most likely related 
to vortex nucleation j|J^-[7|,[l5|,[l6 19 2(J , which is usually 
smaller than the Landau and Feynman critical velocities 
H at which phonons |Lq] or vortices |L7]] become ener- 



getically favorable. 

In conclusion, we have demonstrated a new technique 
to study the onset of dissipation in a stirred condensate. 
This method directly measures density distributions and 
pressure differences. The results obtained are consistent 
with calorimetric measurements but are more immune 
to shot-to-shot noise because multiple non-perturbative 
images of the flow field are taken. This technique is well 
suited to shed further light on contradictory predictions 
about the onset of dissipation in a gaseous superfluid 
p"5| poj and to study in detail how the drag force depends 
on condensate size and density, as well as the object size. 
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